3.1737 \(\int (d+e x)^m (a^2+2 a b x+b^2 x^2)^{5/2} \, dx\)

Optimal. Leaf size=337 \[ -\frac{\sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^5 (d+e x)^{m+1}}{e^6 (m+1) (a+b x)}+\frac{5 b \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^4 (d+e x)^{m+2}}{e^6 (m+2) (a+b x)}-\frac{10 b^2 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^3 (d+e x)^{m+3}}{e^6 (m+3) (a+b x)}+\frac{10 b^3 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^2 (d+e x)^{m+4}}{e^6 (m+4) (a+b x)}-\frac{5 b^4 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e) (d+e x)^{m+5}}{e^6 (m+5) (a+b x)}+\frac{b^5 \sqrt{a^2+2 a b x+b^2 x^2} (d+e x)^{m+6}}{e^6 (m+6) (a+b x)} \]

[Out]

-(((b*d - a*e)^5*(d + e*x)^(1 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(1 + m)*(a + b*x))) + (5*b*(b*d - a*e)^
4*(d + e*x)^(2 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(2 + m)*(a + b*x)) - (10*b^2*(b*d - a*e)^3*(d + e*x)^(
3 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(3 + m)*(a + b*x)) + (10*b^3*(b*d - a*e)^2*(d + e*x)^(4 + m)*Sqrt[a
^2 + 2*a*b*x + b^2*x^2])/(e^6*(4 + m)*(a + b*x)) - (5*b^4*(b*d - a*e)*(d + e*x)^(5 + m)*Sqrt[a^2 + 2*a*b*x + b
^2*x^2])/(e^6*(5 + m)*(a + b*x)) + (b^5*(d + e*x)^(6 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(6 + m)*(a + b*x
))

________________________________________________________________________________________

Rubi [A]  time = 0.141864, antiderivative size = 337, normalized size of antiderivative = 1., number of steps used = 3, number of rules used = 2, integrand size = 28, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.071, Rules used = {646, 43} \[ -\frac{\sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^5 (d+e x)^{m+1}}{e^6 (m+1) (a+b x)}+\frac{5 b \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^4 (d+e x)^{m+2}}{e^6 (m+2) (a+b x)}-\frac{10 b^2 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^3 (d+e x)^{m+3}}{e^6 (m+3) (a+b x)}+\frac{10 b^3 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^2 (d+e x)^{m+4}}{e^6 (m+4) (a+b x)}-\frac{5 b^4 \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e) (d+e x)^{m+5}}{e^6 (m+5) (a+b x)}+\frac{b^5 \sqrt{a^2+2 a b x+b^2 x^2} (d+e x)^{m+6}}{e^6 (m+6) (a+b x)} \]

Antiderivative was successfully verified.

[In]

Int[(d + e*x)^m*(a^2 + 2*a*b*x + b^2*x^2)^(5/2),x]

[Out]

-(((b*d - a*e)^5*(d + e*x)^(1 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(1 + m)*(a + b*x))) + (5*b*(b*d - a*e)^
4*(d + e*x)^(2 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(2 + m)*(a + b*x)) - (10*b^2*(b*d - a*e)^3*(d + e*x)^(
3 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(3 + m)*(a + b*x)) + (10*b^3*(b*d - a*e)^2*(d + e*x)^(4 + m)*Sqrt[a
^2 + 2*a*b*x + b^2*x^2])/(e^6*(4 + m)*(a + b*x)) - (5*b^4*(b*d - a*e)*(d + e*x)^(5 + m)*Sqrt[a^2 + 2*a*b*x + b
^2*x^2])/(e^6*(5 + m)*(a + b*x)) + (b^5*(d + e*x)^(6 + m)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(e^6*(6 + m)*(a + b*x
))

Rule 646

Int[((d_.) + (e_.)*(x_))^(m_)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Dist[(a + b*x + c*x^2)^Fra
cPart[p]/(c^IntPart[p]*(b/2 + c*x)^(2*FracPart[p])), Int[(d + e*x)^m*(b/2 + c*x)^(2*p), x], x] /; FreeQ[{a, b,
c, d, e, m, p}, x] && EqQ[b^2 - 4*a*c, 0] &&  !IntegerQ[p] && NeQ[2*c*d - b*e, 0]

Rule 43

Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d
*x)^n, x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && IGtQ[m, 0] && ( !IntegerQ[n] || (EqQ[c, 0]
&& LeQ[7*m + 4*n + 4, 0]) || LtQ[9*m + 5*(n + 1), 0] || GtQ[m + n + 2, 0])

Rubi steps

\begin{align*} \int (d+e x)^m \left (a^2+2 a b x+b^2 x^2\right )^{5/2} \, dx &=\frac{\sqrt{a^2+2 a b x+b^2 x^2} \int \left (a b+b^2 x\right )^5 (d+e x)^m \, dx}{b^4 \left (a b+b^2 x\right )}\\ &=\frac{\sqrt{a^2+2 a b x+b^2 x^2} \int \left (-\frac{b^5 (b d-a e)^5 (d+e x)^m}{e^5}+\frac{5 b^6 (b d-a e)^4 (d+e x)^{1+m}}{e^5}-\frac{10 b^7 (b d-a e)^3 (d+e x)^{2+m}}{e^5}+\frac{10 b^8 (b d-a e)^2 (d+e x)^{3+m}}{e^5}-\frac{5 b^9 (b d-a e) (d+e x)^{4+m}}{e^5}+\frac{b^{10} (d+e x)^{5+m}}{e^5}\right ) \, dx}{b^4 \left (a b+b^2 x\right )}\\ &=-\frac{(b d-a e)^5 (d+e x)^{1+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (1+m) (a+b x)}+\frac{5 b (b d-a e)^4 (d+e x)^{2+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (2+m) (a+b x)}-\frac{10 b^2 (b d-a e)^3 (d+e x)^{3+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (3+m) (a+b x)}+\frac{10 b^3 (b d-a e)^2 (d+e x)^{4+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (4+m) (a+b x)}-\frac{5 b^4 (b d-a e) (d+e x)^{5+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (5+m) (a+b x)}+\frac{b^5 (d+e x)^{6+m} \sqrt{a^2+2 a b x+b^2 x^2}}{e^6 (6+m) (a+b x)}\\ \end{align*}

Mathematica [A]  time = 0.264481, size = 167, normalized size = 0.5 \[ \frac{\left ((a+b x)^2\right )^{5/2} (d+e x)^{m+1} \left (-\frac{10 b^2 (d+e x)^2 (b d-a e)^3}{m+3}+\frac{10 b^3 (d+e x)^3 (b d-a e)^2}{m+4}-\frac{5 b^4 (d+e x)^4 (b d-a e)}{m+5}+\frac{5 b (d+e x) (b d-a e)^4}{m+2}-\frac{(b d-a e)^5}{m+1}+\frac{b^5 (d+e x)^5}{m+6}\right )}{e^6 (a+b x)^5} \]

Antiderivative was successfully verified.

[In]

Integrate[(d + e*x)^m*(a^2 + 2*a*b*x + b^2*x^2)^(5/2),x]

[Out]

(((a + b*x)^2)^(5/2)*(d + e*x)^(1 + m)*(-((b*d - a*e)^5/(1 + m)) + (5*b*(b*d - a*e)^4*(d + e*x))/(2 + m) - (10
*b^2*(b*d - a*e)^3*(d + e*x)^2)/(3 + m) + (10*b^3*(b*d - a*e)^2*(d + e*x)^3)/(4 + m) - (5*b^4*(b*d - a*e)*(d +
e*x)^4)/(5 + m) + (b^5*(d + e*x)^5)/(6 + m)))/(e^6*(a + b*x)^5)

________________________________________________________________________________________

Maple [B]  time = 0.16, size = 1361, normalized size = 4. \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((e*x+d)^m*(b^2*x^2+2*a*b*x+a^2)^(5/2),x)

[Out]

((b*x+a)^2)^(5/2)*(e*x+d)^(1+m)*(b^5*e^5*m^5*x^5+5*a*b^4*e^5*m^5*x^4+15*b^5*e^5*m^4*x^5+10*a^2*b^3*e^5*m^5*x^3
+80*a*b^4*e^5*m^4*x^4-5*b^5*d*e^4*m^4*x^4+85*b^5*e^5*m^3*x^5+10*a^3*b^2*e^5*m^5*x^2+170*a^2*b^3*e^5*m^4*x^3-20
*a*b^4*d*e^4*m^4*x^3+475*a*b^4*e^5*m^3*x^4-50*b^5*d*e^4*m^3*x^4+225*b^5*e^5*m^2*x^5+5*a^4*b*e^5*m^5*x+180*a^3*
b^2*e^5*m^4*x^2-30*a^2*b^3*d*e^4*m^4*x^2+1070*a^2*b^3*e^5*m^3*x^3-240*a*b^4*d*e^4*m^3*x^3+1300*a*b^4*e^5*m^2*x
^4+20*b^5*d^2*e^3*m^3*x^3-175*b^5*d*e^4*m^2*x^4+274*b^5*e^5*m*x^5+a^5*e^5*m^5+95*a^4*b*e^5*m^4*x-20*a^3*b^2*d*
e^4*m^4*x+1210*a^3*b^2*e^5*m^3*x^2-420*a^2*b^3*d*e^4*m^3*x^2+3070*a^2*b^3*e^5*m^2*x^3+60*a*b^4*d^2*e^3*m^3*x^2
-940*a*b^4*d*e^4*m^2*x^3+1620*a*b^4*e^5*m*x^4+120*b^5*d^2*e^3*m^2*x^3-250*b^5*d*e^4*m*x^4+120*b^5*e^5*x^5+20*a
^5*e^5*m^4-5*a^4*b*d*e^4*m^4+685*a^4*b*e^5*m^3*x-320*a^3*b^2*d*e^4*m^3*x+3720*a^3*b^2*e^5*m^2*x^2+60*a^2*b^3*d
^2*e^3*m^3*x-1950*a^2*b^3*d*e^4*m^2*x^2+3960*a^2*b^3*e^5*m*x^3+540*a*b^4*d^2*e^3*m^2*x^2-1440*a*b^4*d*e^4*m*x^
3+720*a*b^4*e^5*x^4-60*b^5*d^3*e^2*m^2*x^2+220*b^5*d^2*e^3*m*x^3-120*b^5*d*e^4*x^4+155*a^5*e^5*m^3-90*a^4*b*d*
e^4*m^3+2305*a^4*b*e^5*m^2*x+20*a^3*b^2*d^2*e^3*m^3-1780*a^3*b^2*d*e^4*m^2*x+5080*a^3*b^2*e^5*m*x^2+720*a^2*b^
3*d^2*e^3*m^2*x-3360*a^2*b^3*d*e^4*m*x^2+1800*a^2*b^3*e^5*x^3-120*a*b^4*d^3*e^2*m^2*x+1200*a*b^4*d^2*e^3*m*x^2
-720*a*b^4*d*e^4*x^3-180*b^5*d^3*e^2*m*x^2+120*b^5*d^2*e^3*x^3+580*a^5*e^5*m^2-595*a^4*b*d*e^4*m^2+3510*a^4*b*
e^5*m*x+300*a^3*b^2*d^2*e^3*m^2-3880*a^3*b^2*d*e^4*m*x+2400*a^3*b^2*e^5*x^2-60*a^2*b^3*d^3*e^2*m^2+2460*a^2*b^
3*d^2*e^3*m*x-1800*a^2*b^3*d*e^4*x^2-840*a*b^4*d^3*e^2*m*x+720*a*b^4*d^2*e^3*x^2+120*b^5*d^4*e*m*x-120*b^5*d^3
*e^2*x^2+1044*a^5*e^5*m-1710*a^4*b*d*e^4*m+1800*a^4*b*e^5*x+1480*a^3*b^2*d^2*e^3*m-2400*a^3*b^2*d*e^4*x-660*a^
2*b^3*d^3*e^2*m+1800*a^2*b^3*d^2*e^3*x+120*a*b^4*d^4*e*m-720*a*b^4*d^3*e^2*x+120*b^5*d^4*e*x+720*a^5*e^5-1800*
a^4*b*d*e^4+2400*a^3*b^2*d^2*e^3-1800*a^2*b^3*d^3*e^2+720*a*b^4*d^4*e-120*b^5*d^5)/(b*x+a)^5/e^6/(m^6+21*m^5+1
75*m^4+735*m^3+1624*m^2+1764*m+720)

________________________________________________________________________________________

Maxima [B]  time = 1.24332, size = 1069, normalized size = 3.17 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^m*(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="maxima")

[Out]

((m^5 + 15*m^4 + 85*m^3 + 225*m^2 + 274*m + 120)*b^5*e^6*x^6 - 60*(m^2 + 11*m + 30)*a^2*b^3*d^4*e^2 + 20*(m^3
+ 15*m^2 + 74*m + 120)*a^3*b^2*d^3*e^3 - 5*(m^4 + 18*m^3 + 119*m^2 + 342*m + 360)*a^4*b*d^2*e^4 + (m^5 + 20*m^
4 + 155*m^3 + 580*m^2 + 1044*m + 720)*a^5*d*e^5 + 120*a*b^4*d^5*e*(m + 6) - 120*b^5*d^6 + ((m^5 + 10*m^4 + 35*
m^3 + 50*m^2 + 24*m)*b^5*d*e^5 + 5*(m^5 + 16*m^4 + 95*m^3 + 260*m^2 + 324*m + 144)*a*b^4*e^6)*x^5 - 5*((m^4 +
6*m^3 + 11*m^2 + 6*m)*b^5*d^2*e^4 - (m^5 + 12*m^4 + 47*m^3 + 72*m^2 + 36*m)*a*b^4*d*e^5 - 2*(m^5 + 17*m^4 + 10
7*m^3 + 307*m^2 + 396*m + 180)*a^2*b^3*e^6)*x^4 + 10*(2*(m^3 + 3*m^2 + 2*m)*b^5*d^3*e^3 - 2*(m^4 + 9*m^3 + 20*
m^2 + 12*m)*a*b^4*d^2*e^4 + (m^5 + 14*m^4 + 65*m^3 + 112*m^2 + 60*m)*a^2*b^3*d*e^5 + (m^5 + 18*m^4 + 121*m^3 +
372*m^2 + 508*m + 240)*a^3*b^2*e^6)*x^3 - 5*(12*(m^2 + m)*b^5*d^4*e^2 - 12*(m^3 + 7*m^2 + 6*m)*a*b^4*d^3*e^3
+ 6*(m^4 + 12*m^3 + 41*m^2 + 30*m)*a^2*b^3*d^2*e^4 - 2*(m^5 + 16*m^4 + 89*m^3 + 194*m^2 + 120*m)*a^3*b^2*d*e^5
- (m^5 + 19*m^4 + 137*m^3 + 461*m^2 + 702*m + 360)*a^4*b*e^6)*x^2 - (120*(m^2 + 6*m)*a*b^4*d^4*e^2 - 60*(m^3
+ 11*m^2 + 30*m)*a^2*b^3*d^3*e^3 + 20*(m^4 + 15*m^3 + 74*m^2 + 120*m)*a^3*b^2*d^2*e^4 - 5*(m^5 + 18*m^4 + 119*
m^3 + 342*m^2 + 360*m)*a^4*b*d*e^5 - (m^5 + 20*m^4 + 155*m^3 + 580*m^2 + 1044*m + 720)*a^5*e^6 - 120*b^5*d^5*e
*m)*x)*(e*x + d)^m/((m^6 + 21*m^5 + 175*m^4 + 735*m^3 + 1624*m^2 + 1764*m + 720)*e^6)

________________________________________________________________________________________

Fricas [B]  time = 1.84347, size = 3060, normalized size = 9.08 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^m*(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="fricas")

[Out]

(a^5*d*e^5*m^5 - 120*b^5*d^6 + 720*a*b^4*d^5*e - 1800*a^2*b^3*d^4*e^2 + 2400*a^3*b^2*d^3*e^3 - 1800*a^4*b*d^2*
e^4 + 720*a^5*d*e^5 + (b^5*e^6*m^5 + 15*b^5*e^6*m^4 + 85*b^5*e^6*m^3 + 225*b^5*e^6*m^2 + 274*b^5*e^6*m + 120*b
^5*e^6)*x^6 + (720*a*b^4*e^6 + (b^5*d*e^5 + 5*a*b^4*e^6)*m^5 + 10*(b^5*d*e^5 + 8*a*b^4*e^6)*m^4 + 5*(7*b^5*d*e
^5 + 95*a*b^4*e^6)*m^3 + 50*(b^5*d*e^5 + 26*a*b^4*e^6)*m^2 + 12*(2*b^5*d*e^5 + 135*a*b^4*e^6)*m)*x^5 - 5*(a^4*
b*d^2*e^4 - 4*a^5*d*e^5)*m^4 + 5*(360*a^2*b^3*e^6 + (a*b^4*d*e^5 + 2*a^2*b^3*e^6)*m^5 - (b^5*d^2*e^4 - 12*a*b^
4*d*e^5 - 34*a^2*b^3*e^6)*m^4 - (6*b^5*d^2*e^4 - 47*a*b^4*d*e^5 - 214*a^2*b^3*e^6)*m^3 - (11*b^5*d^2*e^4 - 72*
a*b^4*d*e^5 - 614*a^2*b^3*e^6)*m^2 - 6*(b^5*d^2*e^4 - 6*a*b^4*d*e^5 - 132*a^2*b^3*e^6)*m)*x^4 + 5*(4*a^3*b^2*d
^3*e^3 - 18*a^4*b*d^2*e^4 + 31*a^5*d*e^5)*m^3 + 10*(240*a^3*b^2*e^6 + (a^2*b^3*d*e^5 + a^3*b^2*e^6)*m^5 - 2*(a
*b^4*d^2*e^4 - 7*a^2*b^3*d*e^5 - 9*a^3*b^2*e^6)*m^4 + (2*b^5*d^3*e^3 - 18*a*b^4*d^2*e^4 + 65*a^2*b^3*d*e^5 + 1
21*a^3*b^2*e^6)*m^3 + 2*(3*b^5*d^3*e^3 - 20*a*b^4*d^2*e^4 + 56*a^2*b^3*d*e^5 + 186*a^3*b^2*e^6)*m^2 + 4*(b^5*d
^3*e^3 - 6*a*b^4*d^2*e^4 + 15*a^2*b^3*d*e^5 + 127*a^3*b^2*e^6)*m)*x^3 - 5*(12*a^2*b^3*d^4*e^2 - 60*a^3*b^2*d^3
*e^3 + 119*a^4*b*d^2*e^4 - 116*a^5*d*e^5)*m^2 + 5*(360*a^4*b*e^6 + (2*a^3*b^2*d*e^5 + a^4*b*e^6)*m^5 - (6*a^2*
b^3*d^2*e^4 - 32*a^3*b^2*d*e^5 - 19*a^4*b*e^6)*m^4 + (12*a*b^4*d^3*e^3 - 72*a^2*b^3*d^2*e^4 + 178*a^3*b^2*d*e^
5 + 137*a^4*b*e^6)*m^3 - (12*b^5*d^4*e^2 - 84*a*b^4*d^3*e^3 + 246*a^2*b^3*d^2*e^4 - 388*a^3*b^2*d*e^5 - 461*a^
4*b*e^6)*m^2 - 6*(2*b^5*d^4*e^2 - 12*a*b^4*d^3*e^3 + 30*a^2*b^3*d^2*e^4 - 40*a^3*b^2*d*e^5 - 117*a^4*b*e^6)*m)
*x^2 + 2*(60*a*b^4*d^5*e - 330*a^2*b^3*d^4*e^2 + 740*a^3*b^2*d^3*e^3 - 855*a^4*b*d^2*e^4 + 522*a^5*d*e^5)*m +
(720*a^5*e^6 + (5*a^4*b*d*e^5 + a^5*e^6)*m^5 - 10*(2*a^3*b^2*d^2*e^4 - 9*a^4*b*d*e^5 - 2*a^5*e^6)*m^4 + 5*(12*
a^2*b^3*d^3*e^3 - 60*a^3*b^2*d^2*e^4 + 119*a^4*b*d*e^5 + 31*a^5*e^6)*m^3 - 10*(12*a*b^4*d^4*e^2 - 66*a^2*b^3*d
^3*e^3 + 148*a^3*b^2*d^2*e^4 - 171*a^4*b*d*e^5 - 58*a^5*e^6)*m^2 + 12*(10*b^5*d^5*e - 60*a*b^4*d^4*e^2 + 150*a
^2*b^3*d^3*e^3 - 200*a^3*b^2*d^2*e^4 + 150*a^4*b*d*e^5 + 87*a^5*e^6)*m)*x)*(e*x + d)^m/(e^6*m^6 + 21*e^6*m^5 +
175*e^6*m^4 + 735*e^6*m^3 + 1624*e^6*m^2 + 1764*e^6*m + 720*e^6)

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)**m*(b**2*x**2+2*a*b*x+a**2)**(5/2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [B]  time = 1.38734, size = 4316, normalized size = 12.81 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^m*(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="giac")

[Out]

((x*e + d)^m*b^5*m^5*x^6*e^6*sgn(b*x + a) + (x*e + d)^m*b^5*d*m^5*x^5*e^5*sgn(b*x + a) + 5*(x*e + d)^m*a*b^4*m
^5*x^5*e^6*sgn(b*x + a) + 15*(x*e + d)^m*b^5*m^4*x^6*e^6*sgn(b*x + a) + 5*(x*e + d)^m*a*b^4*d*m^5*x^4*e^5*sgn(
b*x + a) + 10*(x*e + d)^m*b^5*d*m^4*x^5*e^5*sgn(b*x + a) - 5*(x*e + d)^m*b^5*d^2*m^4*x^4*e^4*sgn(b*x + a) + 10
*(x*e + d)^m*a^2*b^3*m^5*x^4*e^6*sgn(b*x + a) + 80*(x*e + d)^m*a*b^4*m^4*x^5*e^6*sgn(b*x + a) + 85*(x*e + d)^m
*b^5*m^3*x^6*e^6*sgn(b*x + a) + 10*(x*e + d)^m*a^2*b^3*d*m^5*x^3*e^5*sgn(b*x + a) + 60*(x*e + d)^m*a*b^4*d*m^4
*x^4*e^5*sgn(b*x + a) + 35*(x*e + d)^m*b^5*d*m^3*x^5*e^5*sgn(b*x + a) - 20*(x*e + d)^m*a*b^4*d^2*m^4*x^3*e^4*s
gn(b*x + a) - 30*(x*e + d)^m*b^5*d^2*m^3*x^4*e^4*sgn(b*x + a) + 20*(x*e + d)^m*b^5*d^3*m^3*x^3*e^3*sgn(b*x + a
) + 10*(x*e + d)^m*a^3*b^2*m^5*x^3*e^6*sgn(b*x + a) + 170*(x*e + d)^m*a^2*b^3*m^4*x^4*e^6*sgn(b*x + a) + 475*(
x*e + d)^m*a*b^4*m^3*x^5*e^6*sgn(b*x + a) + 225*(x*e + d)^m*b^5*m^2*x^6*e^6*sgn(b*x + a) + 10*(x*e + d)^m*a^3*
b^2*d*m^5*x^2*e^5*sgn(b*x + a) + 140*(x*e + d)^m*a^2*b^3*d*m^4*x^3*e^5*sgn(b*x + a) + 235*(x*e + d)^m*a*b^4*d*
m^3*x^4*e^5*sgn(b*x + a) + 50*(x*e + d)^m*b^5*d*m^2*x^5*e^5*sgn(b*x + a) - 30*(x*e + d)^m*a^2*b^3*d^2*m^4*x^2*
e^4*sgn(b*x + a) - 180*(x*e + d)^m*a*b^4*d^2*m^3*x^3*e^4*sgn(b*x + a) - 55*(x*e + d)^m*b^5*d^2*m^2*x^4*e^4*sgn
(b*x + a) + 60*(x*e + d)^m*a*b^4*d^3*m^3*x^2*e^3*sgn(b*x + a) + 60*(x*e + d)^m*b^5*d^3*m^2*x^3*e^3*sgn(b*x + a
) - 60*(x*e + d)^m*b^5*d^4*m^2*x^2*e^2*sgn(b*x + a) + 5*(x*e + d)^m*a^4*b*m^5*x^2*e^6*sgn(b*x + a) + 180*(x*e
+ d)^m*a^3*b^2*m^4*x^3*e^6*sgn(b*x + a) + 1070*(x*e + d)^m*a^2*b^3*m^3*x^4*e^6*sgn(b*x + a) + 1300*(x*e + d)^m
*a*b^4*m^2*x^5*e^6*sgn(b*x + a) + 274*(x*e + d)^m*b^5*m*x^6*e^6*sgn(b*x + a) + 5*(x*e + d)^m*a^4*b*d*m^5*x*e^5
*sgn(b*x + a) + 160*(x*e + d)^m*a^3*b^2*d*m^4*x^2*e^5*sgn(b*x + a) + 650*(x*e + d)^m*a^2*b^3*d*m^3*x^3*e^5*sgn
(b*x + a) + 360*(x*e + d)^m*a*b^4*d*m^2*x^4*e^5*sgn(b*x + a) + 24*(x*e + d)^m*b^5*d*m*x^5*e^5*sgn(b*x + a) - 2
0*(x*e + d)^m*a^3*b^2*d^2*m^4*x*e^4*sgn(b*x + a) - 360*(x*e + d)^m*a^2*b^3*d^2*m^3*x^2*e^4*sgn(b*x + a) - 400*
(x*e + d)^m*a*b^4*d^2*m^2*x^3*e^4*sgn(b*x + a) - 30*(x*e + d)^m*b^5*d^2*m*x^4*e^4*sgn(b*x + a) + 60*(x*e + d)^
m*a^2*b^3*d^3*m^3*x*e^3*sgn(b*x + a) + 420*(x*e + d)^m*a*b^4*d^3*m^2*x^2*e^3*sgn(b*x + a) + 40*(x*e + d)^m*b^5
*d^3*m*x^3*e^3*sgn(b*x + a) - 120*(x*e + d)^m*a*b^4*d^4*m^2*x*e^2*sgn(b*x + a) - 60*(x*e + d)^m*b^5*d^4*m*x^2*
e^2*sgn(b*x + a) + 120*(x*e + d)^m*b^5*d^5*m*x*e*sgn(b*x + a) + (x*e + d)^m*a^5*m^5*x*e^6*sgn(b*x + a) + 95*(x
*e + d)^m*a^4*b*m^4*x^2*e^6*sgn(b*x + a) + 1210*(x*e + d)^m*a^3*b^2*m^3*x^3*e^6*sgn(b*x + a) + 3070*(x*e + d)^
m*a^2*b^3*m^2*x^4*e^6*sgn(b*x + a) + 1620*(x*e + d)^m*a*b^4*m*x^5*e^6*sgn(b*x + a) + 120*(x*e + d)^m*b^5*x^6*e
^6*sgn(b*x + a) + (x*e + d)^m*a^5*d*m^5*e^5*sgn(b*x + a) + 90*(x*e + d)^m*a^4*b*d*m^4*x*e^5*sgn(b*x + a) + 890
*(x*e + d)^m*a^3*b^2*d*m^3*x^2*e^5*sgn(b*x + a) + 1120*(x*e + d)^m*a^2*b^3*d*m^2*x^3*e^5*sgn(b*x + a) + 180*(x
*e + d)^m*a*b^4*d*m*x^4*e^5*sgn(b*x + a) - 5*(x*e + d)^m*a^4*b*d^2*m^4*e^4*sgn(b*x + a) - 300*(x*e + d)^m*a^3*
b^2*d^2*m^3*x*e^4*sgn(b*x + a) - 1230*(x*e + d)^m*a^2*b^3*d^2*m^2*x^2*e^4*sgn(b*x + a) - 240*(x*e + d)^m*a*b^4
*d^2*m*x^3*e^4*sgn(b*x + a) + 20*(x*e + d)^m*a^3*b^2*d^3*m^3*e^3*sgn(b*x + a) + 660*(x*e + d)^m*a^2*b^3*d^3*m^
2*x*e^3*sgn(b*x + a) + 360*(x*e + d)^m*a*b^4*d^3*m*x^2*e^3*sgn(b*x + a) - 60*(x*e + d)^m*a^2*b^3*d^4*m^2*e^2*s
gn(b*x + a) - 720*(x*e + d)^m*a*b^4*d^4*m*x*e^2*sgn(b*x + a) + 120*(x*e + d)^m*a*b^4*d^5*m*e*sgn(b*x + a) - 12
0*(x*e + d)^m*b^5*d^6*sgn(b*x + a) + 20*(x*e + d)^m*a^5*m^4*x*e^6*sgn(b*x + a) + 685*(x*e + d)^m*a^4*b*m^3*x^2
*e^6*sgn(b*x + a) + 3720*(x*e + d)^m*a^3*b^2*m^2*x^3*e^6*sgn(b*x + a) + 3960*(x*e + d)^m*a^2*b^3*m*x^4*e^6*sgn
(b*x + a) + 720*(x*e + d)^m*a*b^4*x^5*e^6*sgn(b*x + a) + 20*(x*e + d)^m*a^5*d*m^4*e^5*sgn(b*x + a) + 595*(x*e
+ d)^m*a^4*b*d*m^3*x*e^5*sgn(b*x + a) + 1940*(x*e + d)^m*a^3*b^2*d*m^2*x^2*e^5*sgn(b*x + a) + 600*(x*e + d)^m*
a^2*b^3*d*m*x^3*e^5*sgn(b*x + a) - 90*(x*e + d)^m*a^4*b*d^2*m^3*e^4*sgn(b*x + a) - 1480*(x*e + d)^m*a^3*b^2*d^
2*m^2*x*e^4*sgn(b*x + a) - 900*(x*e + d)^m*a^2*b^3*d^2*m*x^2*e^4*sgn(b*x + a) + 300*(x*e + d)^m*a^3*b^2*d^3*m^
2*e^3*sgn(b*x + a) + 1800*(x*e + d)^m*a^2*b^3*d^3*m*x*e^3*sgn(b*x + a) - 660*(x*e + d)^m*a^2*b^3*d^4*m*e^2*sgn
(b*x + a) + 720*(x*e + d)^m*a*b^4*d^5*e*sgn(b*x + a) + 155*(x*e + d)^m*a^5*m^3*x*e^6*sgn(b*x + a) + 2305*(x*e
+ d)^m*a^4*b*m^2*x^2*e^6*sgn(b*x + a) + 5080*(x*e + d)^m*a^3*b^2*m*x^3*e^6*sgn(b*x + a) + 1800*(x*e + d)^m*a^2
*b^3*x^4*e^6*sgn(b*x + a) + 155*(x*e + d)^m*a^5*d*m^3*e^5*sgn(b*x + a) + 1710*(x*e + d)^m*a^4*b*d*m^2*x*e^5*sg
n(b*x + a) + 1200*(x*e + d)^m*a^3*b^2*d*m*x^2*e^5*sgn(b*x + a) - 595*(x*e + d)^m*a^4*b*d^2*m^2*e^4*sgn(b*x + a
) - 2400*(x*e + d)^m*a^3*b^2*d^2*m*x*e^4*sgn(b*x + a) + 1480*(x*e + d)^m*a^3*b^2*d^3*m*e^3*sgn(b*x + a) - 1800
*(x*e + d)^m*a^2*b^3*d^4*e^2*sgn(b*x + a) + 580*(x*e + d)^m*a^5*m^2*x*e^6*sgn(b*x + a) + 3510*(x*e + d)^m*a^4*
b*m*x^2*e^6*sgn(b*x + a) + 2400*(x*e + d)^m*a^3*b^2*x^3*e^6*sgn(b*x + a) + 580*(x*e + d)^m*a^5*d*m^2*e^5*sgn(b
*x + a) + 1800*(x*e + d)^m*a^4*b*d*m*x*e^5*sgn(b*x + a) - 1710*(x*e + d)^m*a^4*b*d^2*m*e^4*sgn(b*x + a) + 2400
*(x*e + d)^m*a^3*b^2*d^3*e^3*sgn(b*x + a) + 1044*(x*e + d)^m*a^5*m*x*e^6*sgn(b*x + a) + 1800*(x*e + d)^m*a^4*b
*x^2*e^6*sgn(b*x + a) + 1044*(x*e + d)^m*a^5*d*m*e^5*sgn(b*x + a) - 1800*(x*e + d)^m*a^4*b*d^2*e^4*sgn(b*x + a
) + 720*(x*e + d)^m*a^5*x*e^6*sgn(b*x + a) + 720*(x*e + d)^m*a^5*d*e^5*sgn(b*x + a))/(m^6*e^6 + 21*m^5*e^6 + 1
75*m^4*e^6 + 735*m^3*e^6 + 1624*m^2*e^6 + 1764*m*e^6 + 720*e^6)